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Abstract. According to a model of the turbulent boundary layer proposed by the authors, 
in the absence of external turbulence the intermediate region between the viscous sublayer 
and the external flow consists of two sharply separated self-similar structures. The velocity 
distribution in these structures is described by two different scaling laws. The mean velocity u 
in the region adjacent to the viscous sublayer is described by the previously obtained Reynolds- 
number-dependent scaling law <fi = uju„ = Arj a , A = lni?eA + |, a = 2 ln ^ eA , f] = u*y/v. (Here 
w* is the dynamic or friction velocity, y is the distance from the wall, v the kinematic viscosity 
of the fluid, and the Reynolds number ReA is well defined by the data) In the region adjacent to 
the external flow the scaling law is different: <j> — Brf '. The power (3 for zero-pressure-gradient 
boundary layers was found by processing various experimental data and is close (with some 
scatter) to 0.2. 

We show here that for non-zero-pressure-gradient boundary layers, the power j3 is larger 
than 0.2 in the case of adverse pressure gradient and less than 0.2 for favourable pressure 
gradient. Similarity analysis suggests that both the coefficient B and the power (3 depend on 
Re\ and on a new dimensionless parameter P proportional to the pressure gradient. Recent 
experimental data of Perry, Marusic and Jones (l)-(4) were analyzed and the results are in 
agreement with the model we propose. 
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1 Introduction 



The model of the turbulent boundary layer at large Reynolds number proposed by Clauser 
(5) and Coles (6) is widely accepted and used. This model is based on the assumpation that 
the transition from the wall region described by the Karman-Prandtl (7)-(8) universal loga- 
rithmic law to the external flow is smooth. On the basis of our analysis of experimental data 
published over the last 30 years we arrived at a different model (see (9)-(ll)). According to 
our model, if the intensity of turbulence in the external flow is low, then the intermediate 
region between the viscous sublayer and the external flow consists of two self-similar struc- 
tures separated by a sharp boundary. In particular, when the the bilogarithmic coordinates 
lg(/>, lg ?7 are used, where = u/u*, r] — u*y/v, the mean velocity profile in the intermediate 
region has a characteristic form of a broken line ('chevron') (Figure 1). In part I of the 
intermediate region the scaling law for the mean velocity distribution takes the form: 

^ = A (m a . [i] 



In part II one finds a different scaling law: 

u 



■■ B 



(?) • Pi 



The constants A,a,B,j3 can be determined with sufficient accuracy by processing the ex- 
perimental data. According to our model the expressions for A and a are identical to those 
in smooth pipes once the Reynolds number is defined correctly: 

^ = _L lnfleA + |, a=^- [3] 

Here Re^ is an effective Reynolds number for turbulent boundary layer, different from the 
usual Reynolds number Reg based on momentum thickness which is arbitrarily though widely 
used in turbulent boundary layer studies. The test of the validity of our model is the closeness 
of two values of lni?eA, lnitei and lni?e 2 , obtained by solving separately the two equations 
[3] with parameters A and a obtained from experimental data. Differences of less than 2- 
3% were obtained in all cases (see (10), (11) for previous data processing), and therefore we 
proposed to take \a.Re^ as half the sum |(lni?ei + lni?e2). 

In region II the power /3 in the scaling law [2] for the zero-pressure-gradient boundary 
layer was found to be close to 0.2 (with some scatter). In cases of non- zero-pressure-gradient 
boundary layers the values of (3 were found to be significantly different from 0.2. In the 
present Note we perform the similarity analysis for non-zero-pressure-gradient turbulent 
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boundary layer. We find that both the coefficient B and the power (3 depend on Re^ and 
on a new similarity parameter P = ud x p/ pu\. We compare the results of this analysis with 
high quality experimental data by Marusic and Perry(l), Marusic (2), Jones, Marusic and 
Perry (3), Jones (4), and come to the instructive conclusions. 

2 The model and the similarity 
analysis 

According to our model the turbulent boundary layer at large Reynolds numbers consists of 
two separate layers I and II. The structure of the vorticity fields in the two layers is different 
although both are self-similar . In layer I the vortical structure is the one common to all 
developed wall-bounded shear flows and the mean flow velocity is described by relations [1] 
and [3]. In these relations Re\ = UA/v, where A is a characteristic length (12) close to 1.6 
of the height of layer I. 

The influence of viscosity is transmitted to the main body of the flow via streaks sep- 
arating from the viscous sublayer.^] The remaining part of the intermediate region of the 
boundary layer is occupied by layer II where the relation [2] holds. It is well known (see 
in particular instructive photographs in Van Dyke's Album of Fluid Motions (13)) that the 
upper boundary of the boundary layer is covered with statistical regularity by large scale 
'humps' and that the upper layer is influenced by the external flow via the form drag of 
these humps as well as by the shear stress. We have shown in earlier work that the mean 
velocity profile is affected by the intermittency of the turbulence, and as the humps affect 
intermittency it is natural to see two different scaling regions. On the basis of these consid- 
erations we have to determine a set of parameters that determine the coefficient B and the 
power f3 in [2]. One of these parameters must be the effective Reynolds number Re^ which 
determines the flow structure in the layer I and is affected in its turn by the viscous sublayer 
and by layer II. The following parameters should also influence the flow in the upper layer: 
the pressure gradient <9 x p (x is the longitudinal coordinate reckoned along the plate; its origin 
is immaterial), the dynamic (friction) velocity m*, and the fluid's kinematic viscosity v and 
density p. The dimensions of governing parameters are as follows 

ra l M r i L r i L ' r i M mi 
[d xP ] = , [uJ = - , [u] = — , [p] = — . [4] 

1 We note that this mechanism for the molecular viscosity affecting the main body of the flow was proposed 
by L. Prandtl in his discussion of Th. von Karman's lecture (8). It is rather astonishing that this idea was 
never repeated in Prandtl's subsequent publications. 
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The first three have independent dimensions so that only one dimensionless governing pa- 
rameter can be formed: 

P [5] 

put 

Thus the parameters B and (3 should depend on two the parameters Re A and P: 

B = B(Re A ,P) , (3 = (3(Re A ,P) . [6] 



3 Comparison with experimental data 

The data for non-zero-pressure-gradient flows are substantially less numerous than data 
for zero-pressure-gradient flows, and do not allow us yet to construct surfaces B(Re\, P), 
(3(ReA,P)- However the high quality data obtained by Marusic and Perry [(1), recently 
brought to completion via the internet] and Jones, Marusic and Perry [(2), also completed 
on the internet] allowed us to come to some instructive conclusions. The experiments of 
Marusic and Perry (1) were performed for two external flow velocities U: 10 m/s and 30 
m/s. The experiments of Jones, Marusic and Perry (3) were performed for three external flow 
velocities U : 5 m/s, 7.5 m/s, and 10 m/s. The results of the processing of the experimental 
data are presented in Table 1. Here x and Reg are given by the authors of the experiments, 
and A = 2| lnitei - lni?e 2 |/(lni?ei + In .Res). 
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For our subsequent analysis we will use the series corresponding to U = 30 m/s of (2) 
and U = 10 m/s of (4) for the following reasons: in spite of a considerable variation in the 
usual parameter Ree, the effective Reynolds number Re a obtained by the the procedure we 
introduced remains nearly constant and close, for U = 30 m/s (2), to a constant lni?eA = 
10.3, and for U — 10 m/s (4) to a constant lni?eA = 10.2. The mean velocity distribution 
in bilogarithmic coordinates for both series is presented in Figure 2. Thus, we are able to 
obtain, with some approximation, cross-sections of the surfaces B(Re\, P), [3(Re\, P). The 
results corresponding to lnRe\ = 10.3 (adverse pressure gradient) are presented in Figures 
3(a) and 3(b); results corresponding to lni?eA = 10.2 (favourable pressure gradient) are 
presented in Figures 3(c) and 3(d). Note that for large values of the favourable pressure 
gradient we were unable to reveal the second self-similar region. The situation is reminiscent 
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of the disappearance of the second self-similar region in flows with an elevated level of free- 
stream turbulence. We found such a situation previously (10) when we processed the results 
of the remarkable experimental work of P.E. Hancock and P. Bradshaw (14). 

In the papers (l)-(4) the results concerning pressure were presented through a coefficient 

r _ v-Voo 
P ~ \pu* 

where is a constant reference pressure. Therefore we calculated the parameter P using 
the relation d x p = \pU 2 d x C v where the density p cancelled out; the values of all the other 
parameters are available in (2), (4). The values of the parameter P for U = 30 m/s (2) and 
for U — 10 m/s (4) are presented in Table 2. 

Table 2 
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Eliminating the parameter P from relations [6], we obtain: 

B = B(Re A ,(3) . [7] 

This relation is presented in Figure 4 in the form of a dependence of B on ^. We see that 
this dependence is close to linear: 

1 7^ 

B = -f~ 2.80 [8] 

for the data by Marusic (2) {adverse pressure gradient) and 

B = i + 1.43 [9] 

for the data by Jones (4) (favourable pressure gradient) 

For layer I there is also a linear relation between the coefficients A and -, but contrary 
to B = B(jj) this relation is universal. The coefficients in the relation B = B(^) should in 
principle depend on Re\. 
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4 Conclusion 



A new similarity parameter is obtained for the flow in the upper self-similar region of a 
developed non-zero-pressure-gradient turbulent boundary layer. Comparison with experi- 
mental data for nearly constant effective Reynolds numbers revealed simple (close to linear) 
Reynolds number-dependent relations between the parameters of the scaling law for the 
mean velocity distributions in the upper self-similar layer. 

The investigation performed in the present Note and the papers (9)-(12) demonstrated 
that the Reynolds number-dependent scaling law for the velocity distribution across the shear 
flow obtained initially for flows in pipes is valid (with the same values of the constants) for 
the developed turbulent boundary layer flows. This allows us to expect that this scaling law 
reflects a universal property of all developed shear flows. The Reynolds number entering 
the law cannot be selected arbitrarily, for example, as Rcq: it is uniquely determined by the 
flow itself. The simple procedure for the determination of the appropriate Reynolds number, 
which we proposed earlier, has been further validated in the present Note. 

We expect that the same approach will work for more complicated flows: mixing layers, 
jets and wall jets. However, the delicate task of investigating such flows requires high quality 
experimental data which are still lacking. 

The concepts of incomplete similarity and vanishing viscosity asymptotics which we used 
for shear flows lead to plausible results for the local strucxture of developed turbulent flows. 
Here, however, high quality experimental data are very rare, specially for the higher order 
structure functions, where we have conjectured that divergences may occur. 

Acknowledgement. The authors express their gratitude to Professor Ivan Marusic for 
clarification of experimental results. The work was supported by the Applied Mathematics 
subprogram of the U.S. Department of Energy under contract DE-AC03-76-SF00098. 
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Figure Captions 



Figure 1. Schematic representation of the mean velocity profile in developed turbulent 
boundary layer in bilogarithmic coordinates lg lg 

Figure 2. (a) The mean velocity profiles in bilogarithmic coordinates in the series of exper- 
iments of Marusic for U = 30 m/s; adverse pressure gradient . 

(1) Re e = 19, 133 (2) Re e = 16, 584, 

(3) Re e = 14, 208 (4) Re e = 10, 997, 

(5) Reg = 8, 588 (6) Re e = 6, 430. 

The 'chevron' structure of the profiles is clearly seen and regions I and II are clearly 
distinguishable. 

(b) The mean velocity profiles in bilogarithmic coordinates in the series of experiments 
of Jones for U = 10 m/s; favourable pressure gradient . 

(1) Re e = 855 (2) Re e = 1, 122, 

(3) Re = l, 314 (4) Re e = 1, 616, 

(5) Re e = 2, 728 (6) Re e = 3, 032. 

The 'chevron' structure of the profiles is clearly seen for the curves (l)-(4), where 
f3 > a. 

Figure 3. (a) Cross-section of the surface P(Re\ : P), for Re A = 10.3; (2). 

(b) Cross-section of the surface B(Re\,P), for Re A = 10.3; (2). 

(c) Cross-section of the surface (3(Re\,P), for Re\ = 10.2; (4). 

(d) Cross-section of the surface B(Re A) P) ) for Re A = 10.2; (4). 

Figure 4. (a) The dependence -B(^) for Re A = 10.3; 
the straight line corresponds to 1.75//? — 2.8. 

(b) The dependence B(±) for Re A ^ 10.2; 
the straight line corresponds to 1/(3 + 1.43. 
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Figure 1 




Figure 2 b 




Figure 3 a 




Figure 3 b 




Figure 3 c 
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